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ABSTRACT 

The theory and equations for forecasting the wind drift of sea ice, 
as previously developed by several researchers, are being used in the fore- 
casting of that component of sea-ice drift. The existing equations are 
modified and programmed for use with a computer to produce a forecast 
of wind drift, thus permitting rapidly made forecasts over large areas. 
Applicability and limitations, ways of handling other factors, i.e. per- 
manent currents and melting and freezing, and methods for general im- 
provement of the programs are discussed. The programs developed are use- 


ful as research tools and as applied to forecasting. 
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1. Introduction. 

The problem of forecasting the wind drift of sea ice has been im- 
portant to man since he first attempted to navigate the frozen sea. Use 
of the polar sea for supply and commercial routes and as an operating 
area for nuclear submarines makes desirable faster methods of ice fore- 
casting for large areas. Zubov, Shuleikin and Fukutomi [1, 2, 3], among 
others, have investigated the problem of wind drift of sea ice and have 
proposed various formulae for prediction of its movement. The U. S. Naval 
Oceanographic Office (formerly U. S. N. Hydrographic Office) has published 
ice-forecasting techniques, including those for forecasting wind drift, 
which are used by Navy meteorologists and others [4]. 

The Oceanographic Office method employs graphs and tables to obtain 
a forecast by hand calculation. Forecasting for any extensive area is a 
slow process and much subjectivity is involved. Data available for hand 
forecasts are rather sparse, making it difficult to do a complete forecast 
even if time were not a factor. 

A machine method of predicting the wind drift of sea ice would represent 
a significant improvement over present techniques; (1) it would rapidly 
provide a complete wind-drift forecast for an area; (2) continuity of fore- 
casts would result from retention of computed values which may be modified, 
in part, by new observations and carried forward into new forecasts; and 
(3) forecast wind information, to serve as input data, is available via 
numerical methods. 

This paper will present one approach to the numerical forecasting of 
wind drift of sea ice, together with an examination of its applicability 


and limitations. 





Ze Background. 

Many expeditions to the Arctic have contributed to available infor- 
mation on the wind drift of ice, geology, meteorology and oceanography 
of this area; these data have been used by later researchers. In September 
of 1893 Nansen's ship FRAM, built to endure the crushing pressures of ice 
floes, entered the pack-ice in latitude 78-50 N off the New Siberian Islands 
for the express purpose of gathering data on the motion of ice and the currents. 
In August of 1896 after drifting for three years with the ice and reaching 
85-55 N, farthest north for any surface ship, FRAM emerged from the ice 
into the Norwegian Sea east of Greenland completing its gathering of 
scientific data in the Arctic. From 1937 to 1940 the Russian icebreaker 
SEDOV drifted over a route similar to that of FRAM. In the intervening 
years several vessels had completed shorter drifts. 

During 1937 a four-man Russian team under Ivan Papanin drifted for 
nine months on an ice floe, the Papanin Floe; thus, the Russians pioneered 
in placing scientific stations on drifting ice islands [Dyson, 5]. Since 
1950 a number of these drifting stations have been occupied by the United 
States and by Russia. 

Zubov [1] utilized data obtained by the SEDOV eypedition to obtain 
empirical rules for the wind drift of close ice (ice which covers most 
of the surface). His results compared remarkably well with those which 
Nansen derived from data obtained during the drift of FRAM. Zubov 
concluded that over a deep sea, away from the influence of land, the wind- 
caused ice drift is, on the average, at an angle of 29° to the right of the 
wind and at a precdt onic is 1/50 of the speed of the wind. In his treatment 
of the drift of open ice (ice covering between two and six tenths of the 
surface) Zubov shows that the speed of drift is a function of concentration 


and increases with both decreasing concentration and increasing roughness. 
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Shuleikin [6], using dynamic relations and the drift of the Papanin 

Floe compared to winds, derived a theory of ice drift in the open ocean. 

In a later writing [2] he examines his theory further by using data from 

the Papanin and NP-1 floes and obtains constants for his equations. 
Shuleikin's comparatively simple equations consider only ice fields of 10/10 
concentration. 

Since Zubov does not give equations for the speed of ice of varying 
concentration, we refer to Fukutomi [3] for a treatment of wind drift of 
sea ice under these conditions. Fukutomi's formulae are considerably 
more complicated than those of Shule1kin, introducing stresses on the 
vertical surfaces of floes by wind and current and stress of wind on the 
exposed water; Shuleikin considered only stresses on the horizontal ice 
surfaces. This treatment by Fukutomi makes necessary auxillary equations 
and more data on floe dimensions. Further, Fukutomi's constants were 
determined from data obtained in the Sea of Okhotsk and thus differ from 
Shuleikin's constants, based on Arctic Ocean data. However, similarities 
exist in the theories and results. Fukutomi concludes that movement 
predicted by his equations agrees well with Nansen's observations. His 
results support Zubov in that drift speed increases with decreasing con- 
centration. Comparison of Fukutomi's results for 10/10 concentration with 
those of Shuleikin shows rather good agreement, particularly for the angle 
of deviation of ice drift from wind direction. Deviation angle is also 
shown to depend on concentration; this is a reasonable result, since an 
increased wind-induced current in exposed water areas directly affects 
ice movement. 

Review of the above papers and others indicatesthat the several 
theories and methods are mutually consistent; their different results are 


related largely to local differences in ice characteristics (e.g. roughness) 





which are described by certain parameters. Therefore, the method developed 
here uses the equations of Shuleikin modified for concentration by reference 


to Fukutomi. 








3. Models for the wind drift and changes in concentration of sea ice. 


Shuleikin [2] has provided the basic equations required to obtain 


close ice movement. These are: 


tan 


and u = vy 


where: 


El 


¥ = tan B + B/u Gly. 
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Figure 1 depicts the relationships among some of the above quantities. 
In deriving the above equations Sheleikin has assumed an ice concentration 
of 10/10. He further assumed the ice drift is in equilibrium with a constant 
wind, i.e. an unaccelerated motion prevails, and all stress is applied to 
the horizontal surfaces of the ice. 

Fukutomi's equations require a knowledge of: (1) total area of ice, 
(2) area of exposed sea, (3) plane surface area for each floe, (4) vertical 
area exposed to wind pressure for each floe, and (5) vertical area exposed 
to current for each floe. Moreover, one must be able to obtain or estimate 
the energy loss due to eddy currents. Most of these factors are not known 
by observation and are difficult to estimate with any degree of accuracy; 
therefore, it is not practicable to use Fukutomi's equations idrectly. 

The equations are too long to reproduce here, and the reader is referred 
to Fukutomi [3] for further study. 

It is possible, however to approximate one set of Fukutomi's curves 
for wind factor, u/V, versus wind speed, V, with Shuleikin's equations by 
the introduction of a correction factor, E = E + ( as ) C1 + 3V' x .01| 
where c is ice concentration and V' is wind speed 1n knots, (c < 10/10). 
The approximation is shown in figure 2 where the curves according to Fukutomi 
represent complete energy loss due to eddy currents; this corresponds well 
to that condition when area of ice is large compared to ice thickness. The 


resulting equations are: 





tan ¥ = tan B + B/uE (3)2 
Bs oe cos ¥ 
and u = V[(E cos a + cos 8) x 5 coun | E (4). 
wow 


The modification was obtained in a purely arbitrary manner in order to best 


fit Fukutomi's curves over the range of winds normally observed. 
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Now, a method is desired to obtain the rate of change of ice con- 
centration at a point and the resulting concentration at that point for 
some future time. Simpson [7] provides equations, which are readily 
adaptable to numerical methods, for computing these quantities. The 


equation of continuity of ice concentration given by Simpson is: 


Ou Ov. 
OC = Sears oc oc 
a -C ( a EE au ) = ee Ox = a Oy (5) ; 


and the equation for concentration at a given time 1S; 


s Oe 5 


where: 


ac ; . 
ae local change in concentration with time 
c = concentration 


u,= x-component of ice-drift velocity 


i 

v.= y-component of ice-drift velocity 

x, = x-coordinate of location 

Via y-coordinate of location 

am initial time 

rea = time at which new concentration is computed 

o , Oc 

ef = time over which — 1s computed, the time step. 
dt 


Equations (3), (4), (5) and (6) permit one to determine wind drift 
of sea ice and concentration changes resulting from that drift, thus 


completing our models. 





4. Development of the computational programs. 


One of the two programs, ICEFOR1, might be considered as a step to- 
ward the final program developed here, using the models developed in 
section 3. ICEFOR1, reproduced in Appendix A, was written using equations 
(1) and (2) (Shuleikin's equations); it produces wind-drift vectors 
for ice fields with k and Kk. known and uses inputs of latitude and wind 
velocity. While it is not the final form of the wind-drift forecast, it 
does have some utility. 

The next step was to program equations (3) and (4), those modified 
for concentration, to obtain outputs of ice drift similar to the output 
from ICEFOR1. Basically a program segment of ICEFOR1 was modified to 
accept all ice concentrations. Changes in concentration were next con- 


sidered. Equation (5) is put into finite difference form: 





u,-u Vy-V 
where: 
ne = time rate of change of concentration 
c = concentration 
u = x-component of drift velocity 


v = y-component of drift velocity 


D = distance between data points and subscripts indicate data 
points as follows: 


a 
Solution of equation (7) requires values of ice-drift velocity (uo>V5) 


and initial concentration, Coe at each data point for each step. A new 
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concentration, c. + Ac, can then be determined for each data point. One 


Ac 
: at any one point data from 


notes, however, that in order to compute iN 
four other points are required. To satisfy this requirement, and to obtain 
an area forecast, a rectangular grid to approximate an area of the earth's 
surface was devised. The grid is a 546-point array of dimensions 21x26. 
It has an adjustable scale which permits the distance between grid points 
to be changed by use of an appropriate scale factor. Obviously if this 
plane is made too large it no longer approximates an area of the earth's 
spherical surface. However, for an area the size of Baffin Bay, requiring 
a scale factor of 2.0 for representing 30 nautical miles between grid 
points, the approximation to the earth's surface is good and no problems 
arise. The grid permits input of data over an area and further permits 
computation of convergence and advection within that area, as discussed 
in detail later. The program is set up to determine ice-drift vectors, 
changes in concentration and new concentration, point -for-point throughout 
the grid; it is named ICDGRID. 

Several authors, notably Browne and Crary [9] and Fukutomi [3], 
have obtained values for the coefficients of friction Ck, and k in 
equations (2) and (4)) which differ from those used by Shuleikin. Also, 
Wittmann and MacDowell [4] apparently use different coefficients in modi- 
fying Shuleikin's equations fcr Baffin Bay. It is apparent that these 
coefficients do vary, as we might expect, for there are variations in the 
roughness of ice with age (newly-formed, thin ice 1s smoother than old 
ridged ice) and with location. Therefore, the programs developed here 
permit the user to enter values of these coefficients which best fit 
the area of forecast. 


Special problems arise when land boundaries are approached, therefore 


Lt 








the program is written to recognize points of concentration entered as 
2.0 as land points and to ignore these points in computing ice drift. 
Certain of the land points are required, however, in the finite differencing 
process used to obtain dc/dt. In order to treat concentrations near the 
boundaries the program is written to see the concentration of land points 
as 10/10 when the drift is onshore and as 0/10 when the drift has an 
offshore component. In this manner changes in concentration along land 
boundaries are made realistic. 

Program ICEGRID is outlined in Appendix B. It utilizes inputs of 
ice concentration, ice thickness and mean daily surface-wind velocity to 
obtain a five-day wind-drift ice forecast in one-day steps. The outputs 
are grids of ice concentration, direction of movement, and 24-hour distance 


of movement for each of the five days of the forecast. 
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5. Applicability 

The programs developed here are applicable to any geographical area. 

If properly indexed in latitude, ICEGRID will compute over the geographical 
area desired, and ICEFOR] computes for geographical points as directed. 
Accurate forecasting in each case requizes that the coefficients of 
Etiection, Kk. and ko which tend to vary with location, be known. Then, 
given primary inputs of ice concentration, ice thickness and wind velocity, 
proper outputs may be produced. 

There are two basic uses mane programs: forecasting the wind-drift 
component of sea-ice drift and research. Both ICEGRID and ICEFORI1 provide 
forecasts of the wind drift of sea ice, with ICEGRID providing the most 
complete forecast. ICDGRID provides an area forecast and gives, in addition 
to direction and distance of movement, a concentration field: ICEFOR1 produces 
only the ice-drift steering vectors at selected geographical points. 

With regard to research JCEFOR] may be used to obtain paired values 
Or K, and K for an area in which sufficient observational data of actual 
drift are available for comparison with computations. The rapid computa- 
tional ability available in machine methods makes this an extremely advantageous 
tool. Additionally, one may easily extend the possibilities for research 
with both programs to include investigation of permanent currents in 
the frozen sea, to study effects of specific features (i.e. pressure ridges, 
floe size) on the wind drift, and to relate drift to geographical features. 
The fact that ICEGRID is applicable to all concentrations permits its use 
in more general research than ICEFOR1, which is restricted to concentra- 


tions of 10/10. 
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6. Limitations. 

This section discusses the limitations of ICEGRID as a forecasting 
method. Since only a limited portion of the total forecasting problem 
is treated in the theory used here, only that part of the drift directly 
due to the wind is found; components due to permanent currents and concen- 
tration changes due to melting and freezing are ignored. Some of the 
limitations discussed are, therefore, not related to the stated problem, 
that of forecasting the wind drift of sea ice. 

Permanent currents. 

Probably the most important limitation of the program is the failure 
to consider the effects of permanent currents. Wherever such currents 
exist they must affect the total drift of the ice. Therefore, the program 
should be expanded to include these effects. 

One may, of course, subjectively superimpose these effects upon the 
computer output, drawing upon a knowledge of the permanent currents to 
infer the total ice movement. A more direct method of considering these 
currents would be to input permanent current vectors at each grid point, 
adding them vectorially to the ice drift due to wind. The total ice 
drift and the resulting changes in concentration would then be given by 
the computer output. 

Curiously, permanent current data are not always readily available. 
Most current observations are made during periods when little or no ice 
is present, and these would likely be poor approximations for times of 
the year when considerable ice cover is present. Still, in many areas, 
say in the East Greenland Current or in the Labrador Current area, a 
calculation based on any approximation of the permanent currents would 


surely be of more use than one considering wind-drift alone. 
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Melting and freezing. 

A problem perhaps equal to that of permanent currents is that of 
assessing the influence of melting or freezing during the forecast period. 
Assuming no effect may be reasonable early in the melting season; however, 
it is most unrealistic at the height of the melting season and early in 
the freezing season. At these times, ice boundaries may actually move 
opposite to the drift if these effects are large enough. This time of 
year from greatest melting to early freezing is also the mat important 
to the forecaster, as this 1s the time when ships ply the ice covered 
sea and the need for an accurate forecast is greatest. Although the 
problems are complex, models which account for ice melting and freezing 
exist and the necessary meteorological and oceanographic data may be 
extrapolated or approximated from climatology. A machine method of ac- 
cumulating degree days of frost or warming degree days might then lead to 
machine-produced estimates of freezing or melting to be incorporated into 
the ice forecast; thus a more complete forecast would result. 

Hummocking and ridging. 

The effects of additional hummocking and ridging are neglected once 
k, and Kk are chose; actually, the drift characteristics of the ice are 
altered through changes in the roughness parameters. The immediate re- 
sult of failing to consider the above factors is a forecast of slower 
than actual drift in areas characterized by increased roughness; for, 
as previously stated, Zubov [1] concludes that increased hummocking increases 
drift speed. In Fukutomi's development [3] hummocking and ridging are con- 
sidered in that they increase vertical areas exposed to both wind pressure 
and currents. Clearly, these factors should be given some consideration; 
but with presently available data any modification would be largely sub- 


jective, as aerial observations indicate only the presence and not the 
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dimensions of such features; further, one may not know the extent of these 
features present in the data from which the determination of k, and k 
was made. 

It is also highly desirable to be able to predict the formation of 
ridges and hummocks, as ice dominated by these forms is difficult, and per- 
haps impossible, to navigate. Large hummocks and ridges also pose a hazard 
to submarine navigation by restricting water depth. They may lend to in- 
creased formation of fast ice, if they gound in shallow water, and thus 
alter the effective land boundary. This problem requires further study 
directed toward a dynamical solution, but will not be considered here. 

Transport from outside the area. 

No consideration is made of ice advected into the forecast area from 
beyond 1ts boundaries, since no means is provided for entering data from 
outside the area into the machine. This inability to compute new con- 
centrations at the grid boundary points leads to corresponding errors in 
concentration near the boundaries. The problem may be further complicated 
by iterated use of the output as continuity data for new forecasts. If 
boundary points are not modified by new observational data anomalous fore- 
casts near the boundary can result, since the original input concentratidns 
will have been retained at the boundary points. Accordingly, the numerical 
forecaster must be alert to change such points from climatology when such 
changes are required in the interest of continuity and common sense. Of 
course, one might use a grid scale to forecast over an area larger than 
that actually required; thus, the boundary points would have less effect 
on the desired forecast area. 

Nearshore drift. 


In the nearshore region of a water body, including that portion near 
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the edge of the fast ice, the drift of sea ice may be quite unlike that 

in the open ocean. Zubov [1] discusses the nearshore drift with respect 

to the effect of winds in piling up or driving away the ice from the shores. 
Shuleikin compares ice drift along East Greenland with that in the open 
Arctic and shows the differences which exist. From these and other dis- 
cussions, it is apparent that the presence of a solid boundary causes 

local changes in concentration and alters the direction and speed of ice 
Giver & . 

The effect of a solid boundary on concentration is to act as a nucleus 
for increasing concentration when there are onshore winds; and, contrarily, 
to act as a zero concentration point in the case of offshore winds (there 
is no source of sea ice on the land to contribute to offshore advection). 
This portion of the problem is considered in ICEGRID by forcing the computer 
to see concentrations at adjacent land points as 0/10 and 10/10 for offshore 
and onshore winds, respectively, in computing changes in concentration at 
nearshore grid points. 

The influence of the boundaries on the drift vector itself is not con- 
sidered. To make such consideration would -require that the program have 
available not only the location of land points but the orientation of the 
boundary as well. To insert such information into the program would have 
required tying it to a specific geographical location which would have made 
it less general; therefore, no attempt was made to handle this effect. 

Thickness. 

The input of ice thickness is as much a problem as considering hummocks 
quantitatively. Few observations of thickness are available, from an oc- 
Peenona| ship report during the navigable season or from a land station; 


these stations may report regularly on fast ice in the area. The result 
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is that thickness must be input from climatology or long-range ice outlooks 
formed from climatology and heat budget considerations. These provide 
reasonably good data over large areas, but local anomalies may be entirely 
smoothed out, especially by the use of climatology. Fortunately my compu- 
tations using Shuleikin's equations show that drift distance decreases 

on the order of 0.5 nautical miles per day as ice thickness increases from 5 
to 8 feet, for all wind speeds. The angle of deviation from the wind increas- 
es with increasing thickness; however, for a thickness change from 5 to 8 feet 
the maximum difference in the angle is about 15° at low wind speeds and less 
with greater wind speeds, being on the order of 5° for 20-knot winds. Thus, 
the climatology or ice outlook thicknesses will give generally acceptable 
results when used in computations. 

Transient states in the ice field. 

In the equations for sea ice drift it is assumed that the ice moves 
everywhere in equilibrium with the mean wind over the period of a forecast 
step, in this instance one day. It should be apparent that, due to inertia 
of the ice floes, a transient stage occurs whenever a change in wind velocity 
occurs. Zubov, however, states that, "all observers are struck by the speed 
with which ice reacts to any change in the nce indicating transient 
states in the ice movement to be rather short. They are ignored, therefore, 
since the program steps are of one-day duration. It might be necessary to 
consider these states for shorter program steps. 

Observational data. | 

Any forecasting program may be severely limited by the availability, 
accuracy and form of data to be used as input. Currently available obser- 


vations of sea ice limit the forecast in all of these respects. The majority 


azinoee No Ni wepcete tee. Moscow, 1953. pp, 355. 
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of data are from aerial reconaissance flights made at somewhat irregular 
intervals over the limits of the ice; thus the amount of data available 
is adequate over only a small area. 

In addition to data scarcity, data reliability is often open to 
question, since the location of ice observed from the air is directly 
dependat on the accuracy with which the position of the aircraft can be 
fixed; and concentrations depend on the accuracy of the human eye. The 
Leningrad Division of the State Oceanographic Institute has tested an 
instrument designed to measure ice coverage. In the tests, during which 
aircraft position could be accurately determined, it was noted that in- 
dependent visual observers not only differed with observations made with 
the instrument, but among themselves as well. Errors in positioning the 
ice boundary were up to 15 miles; errors in area covered were up to 40% 
[Betin; et al, 11]. The implication is that even direct observations 
may reflect the actual condition inaccurately. 

Current ice reporting messages are not readily adaptable to digital 
methods. A report must be plotted on a chart before data can be extracted 
for program inputs. This is a time-consuming process and results in a 
great loss in the speed advantage realized by a machine method. 

Programming language. 

The program as written is in the programming language called FORTRAN '60. 
This language is simple and easy to use and has speeded the production of 
the program. However, it has some disadvantage in its use; for example, 
all data must be input via punched card to magnetic tape, and computed 
values cannot readily be retained for use in new forecasts. Also, since 
a wind- forecasting program is not available in the language, all wind 


forecasts must be extracted from prognostic charts prior to being input. 
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All this results in loss of time in entering data into the program; the 


entire process requires more than 15 hours for a 600 by 750 nautical-mile 
grid, provided the required ice and weather charts are already plotted. 

Program step size. 

Program step size is only a limitation in that it affects other parts 
of the calculation. We have mentioned briefly transient states of the ice 
field; it is possible that the wind field might be such as to make the 
transient state important even in the 24-hour steps now taken, and these 
stages surely would be important if the steps were shortened to the point 
that steady state did not apply. The step size is important in the com- 
putation of oe, As now programmed this quantity is obtained by a for- 
ward differencing method which can introduce errors over a large number 
of time steps. If the number of time steps were increased, a significant 
error probably would result, necessitating a revision of the method of 
obtaining © . A central-point differencing scheme is suggested as a 


means of accomplishing this if program steps are increased in number. 
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7. Results. 

This section discusses only ICEGRID as ICEFOR1 simply solves Shuleikin's 
equations. Due to the time involved in extracting data for the program input 
only three runs were made for verification. This small number of runs does 
not permit statistical analysis of the forecast versus the observed ice 
field. Analysis is further complicated by the fact that observations often 
cover only a small portion of the area of forecast; hence forecasts were 
made for the few periods when verification data were available. Despite 
these disadvantages, several items of significance are shown by the results. 
The effects of permanent currents, for example, are clearly demonstrated. 
First to be considered, however, is the degree to which we can depend on 
the program as a forecast of the wind-drift component of the sea ice motion. 

In the analysis of results it must be remembered that only the wind- 
induced component of drift is predicted. The immediate conclusion under 
this condition is that ordinarily one would not expect the forecast to 
conform exactly with the observed ice conditions. Such a conformance 
would result only if the drift due to permanent currents and heat budget 
effects were negligible. 

Using coefficients of friction obtained from ICEFOR1, ICEGRID produces 
almost identically the drift which one obtains by hand using the curves 
of the U. S. N. Oceanographic Office [4] for Baffin Bay (these were for 
ag .003 and ian .013). Thus, one is assured that for areas of 10/10 
concentration forecasts of wind drift will be at least as reliable as 
those obtained by hand. Considering that we alter Shuleikin's equations 
for concentration, we may extend this reasoning to say that forecasts 
of wind drift ought to be better than those obtained by hand since hand 


forecasts are made using curves for 10/10 concentration, without modification. 
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While no hand forecast has been made, for such comparison one machine fore- 
cast was done considering all concentrations to be 10/10 for drift calcu- 
lations; a significant drift difference was noted from a forecast modified 
for concentration, using the same input data. The forecast considering 

all concentrations was better for all points in the grid than the one 
assuming ice concentration to be 10/10 everywhere. 

To further establish the worth of ICEGRID in predicting wind drift, 
an examination of the machine forecast versus a forecast based on persis- 
tence of the input field was made. Figures 3, 4 and 5 show the forecast 
boundary versus input and observed boundaries for each of the three fore- 
casts made. In each forecast the machine result was better than persistence, 
in that it gave results which were generally more like the cbserved ice 
field so far as ice boundaries were concerned. The movement was generally 
in the right direction. Comparing forecast concentrations to persistence 
of concentration also shows the forecast to be better since it gives, at 
least, a hint as to the direction of change in concentration. Figures 6 
and 7 contrast an input concentration to its five-day forecast. 

In comparing the machine forecast with the observations at verifying 
time one finds, in general, a considerable discrepancy in the location of 
ice boundaries, say five to 30 miles per five days, and sometimes in con- 
centrations, as much as six tenths, which must arise out of effects other 
than pure wind drift. In most instances, the forecast drift is of the right 
sense, but it is generally insufficient to produce the observed results. 

A full analysis of the forecast results versus observation requires, then, 

a consideration of factors ignored in the development of the method. Most 
of these have been previously discussed (see section 6) and will, therefore, 
be mentioned here only as they seem to affect the validity of a forecast 


made using only wind drift. a 
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The effect of permanent currents is the most apparent of all effects 
exterior to the wind. Forecasts generally show the observed ice boundary 
to vary from that forecast by an increment which, for Baffin Bay, appears 
to conform well with the permanent current believed to exist. It is con- 
ceivable that for extremely light winds the effects of the permenant 
current would dominate leading to a drift entirely different from that 
forecast. In running forecasts, it was noticed that in central Baffin 
Bay around 67 N latitude forecasts conformed rather well with observations, 
indicating a lack of permanent current there, or at least a current parallel 
to the ice boundary. The study of permanent currents through use of this 
program is suggested as a possible research problem. 

A second noticable effect is that due to melting, since the forecasts 
were all made during the melting season. This iS seen most easily through 
comparing observed concentrations with those forecast. Large decreases 
beyond those forecast are sometimes observed. Such changes are likely 
due to melting and not to divergence or advection alone, as no reasonable 
amount of drift could account for them. 

In a previous section data errors were discussed. These contribute 
to apparent errors in the forecast. Since there is, at present, no method 
available to exclude such data errors, one must use all data assuming it 
to be correct. Fortunately, with regard to ice concentration, it is the 
navigability of an area which is important and not the absolute ice con- 
centration. Thus, some error in forecast concentration may be tolorated. 
The forecast concentration has vaJue if it produces an indication as to 
the navigability of the frozen sea. Analysis of program outputs indicate 
the method to be capable of such indication. In the forecasts produced 
it has always been possible to determine rather accurately areas of de- 
creasing concentrations hence, one is able to determine the areas most likely 


to be navigable. 
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Analysis of forecasts indicates the input wind fields to be poorly 
defined with respect to divergence and curl of the wind on the scales 
desired. Thus, small-scale features in the ice field are not shown or are 
smoothed out over the period of forecast. These small features are often 
transient and highly important in the local navigation problem. This raises 
a question as to the effect of small-scale features in the wind field 
on the movement of sea ice. Browne and Crary [9] indicate that the ice 
fields may move with the average wind over a 60-nautical-mile square 
surrounding the location, but are not able to determine conclusively 
this relation, as weather maps fail to supply adequate wind information 
on this scale; and, moreover, all sea ice is not in the form of ice 
fields. It may be, then, that wind forecasts to include smaller-scale 
features in the wind field will be necessary before the smaller-scale 
features of ice drift will become apparent. Since weather reports in 
the Arctic areas are few in number such a definition of the wind field 
may never be possible. 

The results obtained make it apparent that areas of Baffin Bay exist 
wherein nearly all the drift of ice is that due to wind. This seems par- 
ticularly true in the central areanear and just north of the Davis Strait. 
In such an area the wind-drift forecast should constitute a total forecast 
if no heat budget effects are present. 

The few forecasts made using program ICEGRID indicate that the method 
produces good results for wind drift and that it is possible to expand it 


to include other effects. 
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8. Recommendations. 

A review of the preceding sections invites consideration of ways of 
improving the program. One is struck by the need for better methods of 
obtaining and reporting data. Additionally, results obtained point out 
the necessity of expanding the method to provide a complete ice forecast. 

A study of methods of observation and reporting of data has not been 
attempted here and might be the subject of separate investigation. A 
few possible aids to improvement of data quality and availability are 
worthy of note, however. Stereographic aerial photography and TIROS 
(Weather Satallite) photographs are suggested as providing more complete 
observational data. Stereo photographs would give detail as to the dimensions 
of features in the ice such as hummocks, pressure ridges, and even an estimate 
of thickness where ice and open water lie adjacent. Large scale concentra- 
tion data should be available from TIROS photographs where cloud cover does 
not obscure the surface. Scanning of photographs might be done electronically, 
obtaining directly concentration inputs for each grid point. 

Even with better methods of obtaining data, time required to input 
that data to ICEGRID will remain large if the reporting format is not 
changed. The present format presents location, amounts and characteristics 
of the ice verbally, requiring the data to be plotted, extracted for each 
grid point and punched onto data cards prior to computer input. A suitable 
format, from which a direct input to the computer could be obtained, might 
be similar to that of the World Meteorological Organization (WMO) used to 
report weather. Such a report format is suggested in figure 8 and could 
be included in the WMO code. This format would report total concentration, 
concentrations by size and boundary location points for the concentration, 
all reported digitally. Additional plain language data could be included 


at the end of such a message. The code groups could be repeated for each 
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CxCcniipatia. gle lecleat Lol leer 


C @: concentration in tenths 
MN, : tenths of slush, brash, bioke, | 
N, : tenths of small to medium floes 


tals ‘ tenths of large floes & field 

LoL. latitude in whole degrees 

L,L.: longitude in whole degrees as 
| In World Meteorology code 


ic ‘tenths of degrees indicator as 
below: 


Ne jon23| 456 | 789 


FIGURE 8: Suggested ice reporting format 


to be used with the World 
Meteorology code. 
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"concentration area’ of ice observed. The boundary between ice and open 
water could be described by similar latitude-longitude groups following 
a code group, LIMIT. 

We must also seek to provide a means for machine retention of com- 
puted data and a wind-forecasting program. Fortunately, the U. S. Navy 
Fleet Numerical Weather Facility has available a program computing a wind 
prognosis. This suggests rewriting ICEGRID in SCRAP, the programming 
language used by FNWF, Bieaininemeenediate ly a machine forecast of winds 
and, at the same time, the ability to retain computed data as continuity 
for new forecasts. 

Extension of the program to include the effects of permanent currents, 
heat budget effects, land and fast-ice effects on the drift vector (Gn 
short, a complete forecast )is also suggested. Such expansion would require 
considerable time and effort on programming and on obtaining the necessary 
topography and oceanographic data. The production of an efficient and 
useful forecast method dictates the need for this effort, however, and en- 


courages that time be spent on such expansions. 
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9. Conclusions. 

In the foregoing sections the programs presented here have been 
examined and some of the factors foreign to pure wind drift of sea ice, 
which one ought to consider in a complete forecast, have been discussed. 

The programs developed here are concluded to be superior to hand fore- 
casts on speed and on accuracy, due to the modification made for concentra- 
tion and the ability to retain computed values to provide continuity 
to the forecasts. With the inclusion of the effects of permanent currents 
and heat budget, one may expect to obtain a complete and useful ice-drift 
forecast, provided good inputs are available. 

Further, it 1s concluded that the ability to obtain values rapidly for 
the coefficients of friction, Kk, and k leads to the programs' great value 


in research work. 
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APPENDIX A 


PROGRAM ICEFORI1 
Description. 

ICEFOR1 utilizes Shuleikin's equations, without modification, to 
compute the wind drift of sea ice. Provision is made, however, for entry 
of coefficients of friction suitable to the area of computation. Drift 
vectors are computed for one or more locations provided the primary in- 
puts of latitude, ice thickness, wind speed and wind direction are 
entered. Secondary inputs of longitude and date/time group may also be 
entered for reference purposes, but these are not required in any cal- 
culation. 

This program has its primary use as a research tool. Given a set of 
observations of ice drift, corrected for permanent currents, and of corres- 
ponding wind velocities, one may find the Shuleikin curve with the pair 
of coefficients of friction, k, and ko which best fit the data. Having 
obtained the constants, one may obtain ice-drift vectors to be used in 
forecasting. 

Table A-I defines the variables used in the program language. Where 
variables are equivalent to those in Shuleikin's equations symbols are 
included. The complete program is listed following this discussion. 

Prior to use of the program, it is necessary to assign values to the 
constants AKA and WKW in the subprogram, DRIFOR. The values assigned these 
constants in Figure A-1l are for Baffin Bay. The values assigned to RHOAIR 
and DSUBW may also be changed, but this is not usually required. 

Data input is accomplished using IBM cards with punched values of 
variables entered on them. The first card specifies latitude, longitude, 


date/time group, ice thickness and wind velocity, assumed constant for three 


So 





hours. The second and subsequent cards are 
of additional three-hourly wind data at the 
ing the last card specifying wind data is a 
entered as 5000. This signals the computer 


specifying the complete data, i.e.: another 


available for specification 
previous data point. Follow- 
card on which wind speed is 
to proceed to the next card 


"first card". When all data 


have been read, a card which specified latitude as 300 signals the computer 


to read no more data cards. 


Output of results 


The program prints output in columns, titled as follows: 


Result 


latitude 
longitude 
date-time 

ice speed 

ice direction 
distance/3 hrs. 
psi 
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Units 


degrees 

degrees 

day-hour 

knots 

degrees 

nautical miles/day 
radians 





PLAT 
PLONG 
DTG 
THICK 
VELWIND 
DIR WIND 
VELICE 
DIRICE 
DIST3 
AKA 
WKW 
RHOAIR 
DSUBW 
ENN 
OMEGAB 
XT 
ALFA 
BETA 
SQMASS 
BEE 
PSTO 
PSI 
VVICE 


tot wow u @ uu 
wWeDwaut 


TABLE A-lI 
Definitions of variables in ICEFORIL 


latitude 

longitude 

data and time group 

ice thickness 

wind speed in knots 

wind direction 

ice-drift speed in knots 

ice-drift direction 

3-hour drift distance in nautical mile 
coefficient of friction between air and ice, k 
coefficient of friction between ice and water, k 
air density, 6, is 
water density, 6 
N — 

wo sin § =w 


WwW 


an angle, used in computation 
angle between wind and ice-drift direction, ¥ 


= ice-drift speed in meters per second, u 
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APPENDIX B 
PROGRAM ICEGRID 
Description. 

ICEGRID computes the wind drift of sea ice using equations after 
Shuleikin, as modified for concentration less than 10/10, and concentra- 
tion changes according to the equations of Simpson. Drift vectors are 
computed over a 546-point array using latitude, ice thickness, ice con- 
centration, wind speed and wind direction as inputs. The purpose of the 
program is to obtain a five-day forecast of the wind drift of sea ice. 

Table B-I defines the variables used in the program language. Where 
variables are equivalent to those of Shuleikin or Simpson, symbols are 
included. The complete program is listed following this discussion. 

Data input is accomplished using IBM cards with the values of the 
variables punched on them. Four constants must also be entered into the 
program; these are also input via IBM card. 

The first data card specifies the constants, SCALE, which is the 
scale of the grid, and BLAT, the base latitude. A scale of one provides 
grid squares of 60 nautical miles on a side; a scale of two gives grid 
squares 30 nautical miles on a side; and so on. BLAT establishes the 
northern-most boundary of the grid and is always the latitude of that 
boundary plus the inverse of the scale factor. The second data card 
specifies the values assigned to the coefficients of friction, AKA and WKW. 
The fourteen data cards following the two just defined specify the concen- 
tration of ice at all of the 546 grid points, values determined by 
averaging observations over the four adjacent grid squares. Land points 
are entered as 2.0. Following these cards are fourteen cards specifying 


ice thickness over the grid, taken to be the mean climatological thickness 
40 





for the period of forecast. Thickness is entered in feet and is obtained 
in the same manner as concentration for each grid point. Fourteen more 
cards specify wind speed in knots, obtained from prognostic charts, for 
each grid point by averaging over the four adjacent grid squares. Fol- 
lowing this group of cards is a group of 21 cards specifying the dixection 
toward which the wind is blowing. 

The remainder of the input cards consists of groups of 14 and 21 cards 
respectively for wind speed and direction for each of the remaining fore- 
cast days. 


Output of results 


Results are printed in grid form with the outer boundary points of 
the grid omitted. The input concentration grid is reproduced as the first 
output. This is followed by grids of concentration, direction of movement 


and distance for each of the five forecast days. 
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CONCICE 
VELWIND 
DIRWIND 


THICK 


TEMPICE 


VELX 
VELX 


DIST 24 


SCALE 
BLAT 
PLAT 
VELICE 
DIRICE 
DCDT 


RCE, YyYECeE, siCE, XXICE = 


AKA 
WKW 
RHOAIR 
DSUBW 
OMEGAB 
XI 

ENN 
ALFA 
BETA 
SQMASS 
BEE 
PSIO 
PSI 
CONFAC 
VVICE 
C 


Ui 


I 


I 


<= 
—- 


= coefficient of friction between ice and water, k 


TABLE B-I 
Definitions of variables in ICEGRID 


ice concentration, c 

wind speed in knots, V' 

wind direction +180 , in degrees 

ice thickness, enters into computation of mass 

a temporary storage for new concentration 
x-component of ice speed, u, in knots 

y-component of ice speed, r, in knots 

distance of ice drift in 24 hours, in nautical miles 
scale factor, S,determines grid size 

base latitude permits proper orientation of grid in latitude 
latitude,degrees 9g 

ice drift speed in knots 

ice drift direction, in degrees 

change in concentration with time, dc/dt 

temporary concentration storages 
coefficient of friction between air and water, k 


W 
air density, 5, 


= water density, 6 


touu tt @ 
WE DwWR Au 


w sin $, BD 


an angle , used for computation 

angle between wind and ice direction, } 

factor modifying for concentration 

ice drift speed in meters per second, u 

constant used to modify outputs of ¥ and u, if consistently 
low or high, normally set to one 
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